Ultrafast dynamical path for the switching of a ferrimagnet after femtosecond heating 
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Ultrafast laser-induced magnetic switching in rare earth, transition metal ferrimagnetic alloys 
has recently been reported to occur by ultrafast heating alone. Using atomistic simulations and 
a ferrimagnetic Landau-Lifshitz-Bloch formalism, we demonstrate that for switching to occur it is 
necessary that angular momentum is transferred from the longitudinal to transverse magnetiza- 
tion components. This dynamical path leads to magnetization switching and subsequent ultrafast 
precession caused by the inter-sublattice exchange field on the nanoscale. 



The behavior of magnetization dynamics triggered by 
an ultrafast laser stimulus is a topic of intense research 
interest in both fundamental and applied magnetism [JJ. 
A range of studies using ultrafast laser pulses have shown 
very different timescales of demagnetization for different 
material; from 100 fs in Ni to 100 ps in Gd [3 ]. For fu- 
ture applications one should aim not only at ultrafast 
demagnetization, but also to controlled magnetization 
switching. 

Magnetization reversal induced by an ultrafast laser 
pulse has been reported in the ferrimagnet GdFeCo, to- 
gether with a rich variety of phenomena [3HE]- Several 
hypotheses have been put forward to explain the ob- 
served magnetization switching: crossing of the angular 
momentum compensation point @], the Inverse Faraday 
Effect j5] , and its combination with ultrafast heating [6] . 
It has been shown that the rare earth (RE) responds 
more slowly to the laser pulse than the transition metal 
(TM) [7], even though the sublattices are strongly ex- 
change coupled. Intriguingly, Radu et. al. [7] show ex- 
perimentally and theoretically the existence of a transient 
ferromagnetic-like state which exists for a few hundred 
femtoseconds. Recently [5], the atomistic model outlined 
in [7J |S] predicted the phenomenon of magnetization re- 
versal induced by heat alone, in the absence of any ex- 
ternal field; a prediction verified experimentally. This 
remarkable result opens many interesting possibilities in 
terms of ultrafast magnetization reversal and potential 
areas of practical exploitation, however a complete theo- 
retical understanding of this effect is currently missing. 

In magnets consisting of more than one magnetic 
species, excitation of the spins on a time scale compa- 
rable with that of the inter-sublattice exchange takes the 
sublattices out of equilibrium with each other. It is in 
this regime where the thermally driven switching of fer- 
rimagnetic GdFeCo occurs. A recent study by Mentink 
et al. [TU] proposed an explanation of the process us- 
ing a model of the magnetization dynamics, a feature of 
which is the separation of the relaxation mechanisms into 
two timescales: one governed by the inter-sublattice ex- 
change and another by the relativistic contribution (cou- 
pling to external degrees of freedom). However, in ultra- 
fast processes the relativistic contribution is most proba- 



bly dominated by electron-electron scattering [TJJ and/or 
by enhanced spin-orbit interaction |12j. both acting on 
the timescale of the exchange interactions. Consequently, 
the two contributions (exchange and intrinsic damping) 
cannot a-priori be separated. This is also confirmed by 
previous atomistic simulations |13j. where it was shown 
that on the fs-ps timescale the longitudinal relaxation of 
ferromagnets is simultaneously defined by the coupling 
to the external bath and the exchange interaction. Ad- 
ditionally, the model in ([TO]) uses the Onsager principle, 
generally valid in equilibrium only. The model is based 
on the physically plausible argument that the switching 
is driven by angular momentum transfer due to the ex- 
change between sub-lattices and the assumption of a lin- 
ear path to reversal, thereby neglecting the role of any 
perpendicular component. Thus far, a complete explana- 
tion of the heat driven, ultrafast reversal process remains 
elusive. 

In this Letter we present a general formalism, leading 
to a macroscopic dynamical equation for a ferrimagnet. 
This is in the form of a Landau-Lifshitz-Bloch (LLB) 
equation, in which the two relaxation mechanisms are 
parts of the same relaxation term. Supported by atom- 
istic calculations, we show that the role of the perpen- 
dicular magnetization component in the ultrafast, heat- 
induced magnetization reversal process is essential, and 
that without it reversal cannot occur. Predictions of the 
ferrimagnetic LLB equation are supported by calcula- 
tions using an atomistic model based on the Heisenberg 
exchange Hamiltonian with Langevin Dynamics. 

In the absence of any external stimulus, the energet- 
ics of the atomistic spin model are described purely by 
exchange interactions, given by the spin Hamiltonian: 

H = - £ • / v H - • S - C 1 ) 

where J™ is the exchange integral between spins i and 
j (i,j are lattice sites), Sj is the normalized magnetic 
moment |Sj| = 1. We model the magnetization dynamics 
of the system using the Landau-Lifshitz-Gilbert (LLG) 
equation with Langevin dynamics, as detailed in Ref. HI 
To simulate the effect of an ultrafast heat pulse we use 
a step-like temperature pulse of duration 500fs with a 
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value of T = T max . In order not to mix different effects, 
we assume the same electron temperature and the same 
damping constants for Gd and Fe spin moments. The 
slow character of the demagnetization in Gd will result 
from the large moment of Gd relative to Fe. Atomistic 
models have proven to be a powerful tool in predicting 
heat-induced switching [S], but fail to provide a simple 
picture for the cause. 

The magnetization dynamics of ferrimagnets is of- 
ten described by two coupled temperature-dependent 
macrospins within LLG micromagnetic equations. How- 
ever, even for a ferromagnet, it has been shown that the 
LLG approach is not a valid micromagnetic equation for 
high temperatures and it has been superceded by the 
LLB equation |13j . Up until now, the LLB equation ex- 
isted only for a single species ferromagnet. In this letter 
we present the LLB equation for a two species system. 
The detailed derivation [H] makes use of the mean-field 
approach (MFA) based on the Heisenberg Hamiltonian 
Eq. [I] and is similar to the original derivation of the LLB 
equation for a ferromagnet [T5]. Importantly, unlike the 
approach used in RefflO| the derivation does not use the 
Onsager principle and is thus valid far from equilibrium. 
Our equation also does not separate the timescales of 
different relaxation sources. 

The resulting ferrimagnetic LLB equation describes 
the average magnetization dynamics in each sublattice 
m„ = ( s i), where v stands for TM or RE sublattice and 
i for spins in the sublattice v . In the absence of an ap- 
plied field, the LLB equation for the TM is written as: 
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with a complementary equation for the RE. The ex- 
change field from the RE is calculated via the MFA [5] 
as Hg X = (zx\ Jtr|/^t)hir, where Jtr is the inter- 
sublattice exchange parameter, z is the number of near- 
est TM neighbors in the ordered lattice, x is the RE 
concentration, fix is the TM magnetic moment, 7t is 
the gyromagnetic ratio for TM, tr = |(mx • iiir)|/tot 
is the absolute value of the projection of the RE mag- 
netization onto the TM magnetization and r e ^ is its 
equilibrium value, a T (T) = 2At&bT'/Jo i t and asp(T) = 
At(1 — fer/J 0j T) are the temperature-dependent TM 
longitudinal and transverse damping parameters and At 
is intrinsic coupling to the bath parameter. Finally, 
Jq.t = (^o,T"i e ,T - ■/o,TR"»e,R)K,T, where m e>T and 
m ej R, are the equilibrium magnetization of the TM and 
RE, Jo,t = (1 — x)z Jt and Jo,tr = xzJth are the zero 
Fourier components of the TM intra and inter-sublatticc 
exchange interaction [9]. Note that even if the micro- 
scopic damping of the TM and RE are the same, the 
macroscopic damping is different due to a different tem- 
perature dependence of magnetization. The third and 



fourth terms in Eq.([2| define the longitudinal relaxation, 
comprised of; the difference between the longitudinal re- 
laxation of the TM to its equilibrium value, and to the 
equilibrium value of the RE respectively. 

The rate parameters Ttt = A^ T , V 
temperature-dependent and have the following forms 
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where we have used the longitudinal susceptibilities 
XT, || = (dm T /dH) H ^ , Xr,\\ = (dm R /dH) H ^ in inter- 
sublattice and intra-sublattice exchange. In Fig. [I] we 
present the temperature dependence of relaxation rates 
evaluated for the parameters of GdFeCo. The temper- 
ature dependence of the relaxation rates are significant 
for ultrafast switching, since close to T c the sign of the 
longitudinal relaxation term (the third term in Eq. 2) 
changes, such that during the ultrafast heating process 
dm T /dt < goes to dm T /dt > 0. However, and unlike 
the statement in Ref. [10] . the longitudinal relaxation it- 
self cannot change the direction of the vector hit, due 
to the multiplication of the longitudinal relaxation term 
in Eq. ^ by hit- To illustrate this, we consider the re- 
duction of Eq. ([2]) to linear reversal processes only, by 
removing the precessional term. Irrespective of the sign 
of dm T /dt its value must go to zero when the sublat- 
tice magnetizations are individually in equilibrium. Since 
m ej T can never be negative without reversal of the RE, 
no reversal of the TM prior to reversal of the RE is pos- 
sible by this route. The change of sign of dm T / dt results 
because the polarization of the TM by the RE leads to 
a small but positive value of mT- Evidently, although 
the longitudinal magnetization process contributes to the 
timescale of reversal it does not drive the switching pro- 
cess. In order to understand the switching mechanism we 
therefore need to consider both longitudinal and trans- 
verse relaxation. 

In the following we reduce the LLB equation ^ to 
a dynamical system, based on information from atom- 
istic modelling. The phase portrait of this dynamical 
system is used to demonstrate the role of transverse 
fluctuations in the heat-driven reversal process. We 
know from numerical simulations that close to the re- 
versal of the sublattices the RE magnetization mP R is 
larger than that of the TM, so that we can assume that 
, m?r <C m a R where m° D is the abso- 
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lute value of the z-component of the RE magnetization 
at the instant of reversal. In this approximation the lon- 
gitudinal field is positive h\\ > 0: hj ~ ajj^- for the 
case before the heat pulse is removed (tot > rn e ^T = 0) 
and hj ~ (a T /2) [IVr — Ttr] after the heat pulse (tot < 
Tn e Tj T = K). The LLB equation for the TM is reduced 
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FIG. 1. Longitudinal relaxation rates as a function of temper- 
ature in the LLB equation, evaluated for GdFeCo parameters. 
The dashed line shows the TM-RE relaxation rate and the 
solid line is that of the TM-TM interaction. At low tempera- 
tures Ttt ^> Ttr, due to the small value of the susceptibility 
XT, ||) therefore the relaxation of the TM magnetization is al- 
ways to its own equilibrium. However, at temperatures close 
to T c , Ttt < Ttr, thus the TM prefers to relax towards the 
RE magnetization in this regime. 



to the following dynamical system: 
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where p — (tj^) 2 = (m^) 2 -I- (m^,) 2 is the TM transverse 
magnetization component and Qt(T) — zxrrv^ Jtr|/Vt 
is the precessional frequency of the anti-ferromagnetic 
exchange mode. 

In Fig. § Eqs. @ are solved as a function of time 
and are presented as a phase space plot with comparison 
to atomistic simulations. The trajectory = cor- 
responds to a linear dynamical mode and according to 
Eq.fptl) is trivially stable for hj < (normal relaxation 
describing the linear magnetization decrease) and unsta- 
ble in the opposite case (hj > 0). The linear path does 
not lead to reversal; a small initial perpendicular compo- 
nent is required to reverse the direction of the equilibrium 
magnetization. A standard analysis of the dynamical sys- 
tem jij with hj > shows that for ^ and close to 
« 0: dm^/dt < 0. This corresponds to the reversal 
trajectories shown in Fig. [2j The physical interpretation 
is that in this case very small perturbations from = 
will lead to the development of a perpendicular magne- 
tization component, as is indeed observed by numerical 
integration of the dynamical system Q. The correspond- 
ing dynamical trajectories of the TM cross = lead- 
ing to the change of the magnetization direction of the 
TM. The dynamical path indicates the increase of m^, 
i.e. the further development of the precessional motion. 

In the same approximation, the LLB equation for the 
RE reads: 
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FIG. 2. a. 

different initial conditions very close to p 
ditions with very small components of are supposed, b. 
Real trajectories calculated via atomistic simulations for TM 
and RE at different pulse T max temperatures. Small circles 
indicate the point where the heat pulse is removed. 



where the upper sign corresponds to the equation for 
and the lower sign for the one, Or = z(l — 
x)m^\ JtrI/Vr an d h$ > is the RE longitudinal field. 
Thus, the perpendicular motion of the TM triggers the 
corresponding precessional motion of the RE via the an- 
gular momentum transfer (the first two terms of Eq. |2| , 
i.e. via perpendicular components) with the same fre- 
quency f^x, but different amplitude. However, the lon- 
gitudinal change in the RE is slower than that of the 
TM. During this dynamical process in some time in- 
terval the RE and TM have both the same sign of the 
m z -component, forming the transient ferromagnetic-like 
state seen experimentally [7J. Note that the subsequent 
precession has a frequency which is proportional to the 
exchange field and thus is extremely fast. The motion of 
the TM around RE direction and vice versa occurs dur- 
ing and after the ferromagnetic-like state until the system 
has relaxed to equilibrium. 

Following the analysis of the LLB equation, one should 
note that in order for magnetization switching to occur, 
a small 'kick' from the unstable linear motion, producing 
an initial torque, is essential. Thus the direct numeri- 
cal integration of Eq. ^ with initial antiparallel align- 
ment of the RE and TM follows the linear magnetization 
recovery path and does not produce switching. Con- 
versely, the presence of a small angle between the two 
sub-lattices provides the necessary condition for magne- 
tization switching, as demonstrated in Fig. [3] It is impor- 
tant to emphasize that in reality non-collinearity of the 
sublattices is always present because of the high temper- 
atures associated with the reversal mechanism. This is 
true in the full atomistic approach with stochastic fields 
as well as the integration of the stochastic LLB equa- 
tion [IB] . 

An outstanding question is whether the magnetiza- 
tion precession, a central prediction of our model, can 
be observed experimentally on a macroscopic sample. 
We should recall that at high temperatures the corre- 
lation between atomic sites is weak, thus we cannot ex- 
pect the precession to occur with the same phase in the 
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FIG. 3. Numerical integration of the switching behaviour for 
the non-stochastic LLB with a small angle between sublattices 
(solid lines). Without the angle (dashed lines) switching does 
not occur, as predicted. The time t = corresponds to the 
end of the laser pulse. 



whole sample; an effect which would make the preces- 
sion macroscopically unobservable. To demonstrate the 
effect we present in Fig. [4] the results of atomistic switch- 
ing simulations in GdFeCo for different system sizes. We 
clearly observe that for small system sizes transverse os- 
cillations with the frequency of an exchange mode are 
visible, consistent with the prediction of our analytical 
model. However, in large system sizes of the order of 
(20 nm) 3 it is averaged out, consistent with the exci- 
tation of localized exchange modes with random phase. 
Note that the same effect happens for very high tem- 
peratures where the observed magnetization trajectory 
appears close to linear; although we stress again the im- 
portance of a small perpendicular component to initiate 
the magnetization reversal, which will occur on a local 
level as demonstrated by Fig. [4] 

The above analysis also allows us to have a clearer 
understanding of the ferromagnetic-like state. The ve- 
locity of the motion is proportional to a^Q,? ~ Jtr- 
Thus, to create a longer transient ferromagnetic state, 
we should decrease the precessional frequency, which can 
be achieved by decreasing the inter-sublattice exchange 
Jtr- In addition, the dynamical response of the RE can 
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FIG. 4. System size dependence of the transverse magnetiza- 
tion components of the TM under ultrafast switching, showing 
cancelation of the localized transverse magnetization compo- 
nents arising from exchange precession for larger system sizes. 
The time t = corresponds to the end of the laser pulse. 
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FIG. 5. Atomistic simulated dynamics of the TM (red dashed 
line) and RE (blue solid line) under the action of an ultrafast 
heat pulse where T =2000 K for 600 fs duration, which is 
then reduced back to OK. The system size is (6 nm) 3 . The 
moment t = corresponds to the end of the laser pulse. The 
inter-sublattice parameter Jtr and Xre were decreased to - 
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J and 0.001 respectively, a. The TM reverses 



very quickly setting up the ferromagnetic-like state which ex- 
ists for around 47 ps. b. Coupled, large amplitude precession 
of the sublattices is clearly seen. 



be tuned by altering the coupling parameter, and so for 
slower dynamics we set \re = 0.001, consistent with a 
small coupling parameter of Gd |17| . Fig. [5] shows sim- 
ulations of the resulting dynamics, and demonstrate an 
extended ferromagnetic state. 

In conclusion, we have outlined a MFA treatment of the 
multispin Heiscnbcrg model as a function of temperature 
resulting in a LLB equation for a ferrimagnet. This equa- 
tion describes the mutual relaxation of sublattices which 
occurs simultaneously under internal damping and inter- 
sublattice exchange. This model allows us to present a 
simple picture of the magnetization reversal of GdFeCo 
in response to an ultrafast heat pulse alone. The physical 
origin of this effect is revealed within the LLB equation 
as a dynamical reversal path resulting from the instabil- 
ity of the linear motion. To trigger the reversal path a 
small perpendicular component is necessary. In practice 
this will arise from random fluctuations of the magnetiza- 
tion at elevated temperatures. The perpendicular com- 
ponent grows in time resulting in ultrafast magnetization 
precession in the inter-sublattice exchange field, also ob- 
served in atomistic simulations for small system sizes. 
The switching is initiated by the TM which arrives at zero 
magnetization faster than the RE and responds dynam- 
ically to its exchange field. Thus, the non-equivalence of 
the two sub-lattices is an essential part of the process. 
Switching into the transient ferromagnetic state occurs 
due to large-amplitude precessional motion of the TM in 
the exchange field from the RE and the slow dynamics of 
RE. The control of the transient ferromagnetic-like state 
can be achieved through tuning of the damping and inter- 
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sublattice exchange in agreement with the dynamics of 
the LLB equation. 
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